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Abstract 

It is well known that a suggestive relation exists that links Schrodinger's equation (SE) to the 
information-optimizing principle based on Fisher's information measure (FIM). We explore here 
an approach that will allow one to infer the optimal FIM compatible with a given amount of prior 
information without explicitly solving first the associated SE. This technique is based on the virial 
theorem and it provides analytic solutions for the physically relevant FIM, that which is minimal 
subject to the constraints posed by the prior information. 



KEYWORDS: Information Theory, Fisher's Information measure, Legendre transform, Virial the- 
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I. INTRODUCTION 



Although Fisher's information measure (FIM) / dates from the 20's, it has attracted intense 



rom physicists only since the mid-90's 



l|-|3|. Interest in FIM's physical appli- 



has been growing exponentially since the appearance of Frieden and So: 



10 



;er s 



24|. 



attention 
cations |3| 

seminal paper jsl- A very small, and certainly not exhaustive sample is that of Refs. 
The mathematical problem of extremizing / under given constraints has proved to be relevant 
in connection with several scenarios related to quantum mechanics and statistical physics, 
as the above cited references indicate. We intend to exhibit novel links between I and the 
Schrodinger equation that, via the virial theorem, lead in natural fashion to a differential 
equation for /. Such equation, that has an analytical solution, encodes the available prior 
knowledge concerning the system at hand in terms of adequately selected expectation values. 
Our solution adds to the rather large Fisher literature a general, explicit expression for that 
particular FIM /Mm that arises out of any constrained J— extremization problem. 

To better understand why this is of importance one should recall that Fisher's information 
and Shannon's entropy play complementary roles 4]- The former is convex, the later concave. 
When one grows, the other diminishes, etc. The associated Shannon's MaxEnt problem has 
as its solution, always, an exponential form that contains those physical quantities whose 
mean values are a priori known. The solution of the FIM minimization problem is instead a 
Schrodinger hke differential equation^, 3,0], whose solutions exhibit a panoply of different 
mathematical forms. We provide here special definite forms for iMm-, expressed in terms of 
those mean values that are a priori known, filling thus a gap in the literature of the physics 
of information. We review basic materials below and start with our presentation in Section 
III-IV. We give illustrative examples in Section V and discuss particular issues in Section 
VI. Our conclusions are given in Section VII. 



II. REVIEW OF BASIC IDEAS 

We briefly review here the formalism developed in Ref. 6] (see also Consider a system 
that is specified by a physical parameter 6 and let /(x, 9) describe the normalized probability 
distribution function (PDF) for this parameter. If an observer were to make a measurement 
of X and had to best infer 9 from such measurement, calling the resulting estimate 9 = 9{x), 
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one might well wonder how well 6 could be determined. Estimation theory asserts that 
the best possible estimator 6 (x), after a very large number of x-samples is examined, suffers a 
mean-square error from 6 obeying the rule le^ = 1, where the Fisher information measure 
(FIM) /, a functional of the PDF, reads 

/ = / dxf{x,e)!^^ ln[/(x,^)]| . (1) 

Any other estimator must have a larger mean-square error (all estimators must be unbiased, 
i.e., satisfy (^(x)) = 6). Thus, FIM has a lower bound. No matter what the parameter ^ of 
the system might be, / has to obey 



Iei>l, (2) 



the celebrated Cramer-Rao bound j^. The particular instance of translational families 
merits a word. They are mono-parametric distribution families of the form f{x, 6) = f{x—6), 
known up to the shift parameter 6. All family members exhibit identical shape. After 
introducing the amplitudes ip such that f{x) = ip{xY, FIM adopts the simpler aspect 



d 
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rfx /(x) |— ln[/(x)]| =Ajdx[ip\x)Y; {dij/dx = ij'). (3) 

Note that for the uniform distribution f{x) = constant one has / = 0. Focus attention now 
a system that is specified by a set of M physical parameters /i^. We can write fik = {^k) with 
Ak = Ak{x). The set of /i^- values constitutes the prior knowledge. It represents available 
empirical information. Let the pertinent probability distribution function (PDF) be f{x). 
Then, 

{Ak) = J dx Ak{x) fix), fc = l,...,M. (4) 

In this context it can be shown (see for example that the physically relevant PDF 

f{x) minimizes the FIM (j3]) subject to the prior conditions and the normalization condi- 
tion. In the celebrated MaxEnt approach of Jaynes' [25| one maximizes the entropy, that 
behaves information- wise in opposite fashion to that of Fisher's measure js]. Normalization 
entails / dxf{x) = 1, and, consequently, our Fisher-based extremization problem adopts the 
appearance 

^{l-^j /(^) - E ^kj dx Ak{x) /(x)j = (5) 

where we have introduced the (M -|- 1) Lagrange multipliers (Aq = a). In Ref. {g] on can 
find the details that lead from (|5]) to a Schrodinger's equation (SE) that yields the desired 
PDF in terms of an amplitude tplx) referred to above before Eq. ([3]). This SE is of the form 



1 9^ 



M 



2 dx- 



fc=i 



(6) 



which can be formally interpreted as the (real) Schrodinger equation for a particle of unit 
mass {% = 1) moving in the effective, "information-related pseudo-potential" ^| 



M 



u 



U{x) 



(7) 



k=l 



in which the normalization-Lagrange multiplier (a/8) plays the role of an energy eigenvalue. 
The Afc are fixed, of course, by recourse to the available prior information. Note that il){x) 
is always real in the case of one- dimensional scenarios, or for the ground state of a real 



potential in N dimensions 



26|. In terms of the amplitudes ipi^x) we have 



jdxf 



'd\nf 

dx 



dx tpl 



dx 



4 / i^n-Q^'ipn dx 



dx 

M 



dx 



, dx^ 



J 'ipn(^a + Afc ipn dx, 



I.e., 



M 

1= a + J2 \k {Ak) . 
fe=i 



(8) 



a form that we will employ in our developments below. The connection between our vari 



ational solutions / and thermodynamics was established in Refs. 6| and 



i i 



in the guise 



of reciprocity relations that express the Legendre-transform structure of thermodynamics. 
They constitute its essential formal ingredient and were re-derived a la Fisher in jo] by 
recasting ([8]) in a fashion that emphasizes the relevant independent variables 



M 



I {{A,),..., {Am)) = a + E {Ak) ■ 

k=l 

The Legendre transform changes the relevant variables. In the case of I this is 



(9) 



M 



a 



I{{Ai) , . . . , {Am)) -Y. {Ak) = a(Ai, . . . , Am), 



(10) 



k=l 



SO that we encounter the three reciprocity relations proved in 
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the last one being a generalized Fisher-Euler theorem. 



III. PRESENT CORE-RESULTS 



Our previous Fisher considerations lead to a scenario for which 



1 (9^ 1 n 

^ = -2^+f^(^)> f/(x) = --EA.A, , E„ = -. (12) 

Enters here, as essential new ingredient in the present considerations, the celebrated virial 

29]. This theorem is inti- 



theorem [28[ that of course applies in this Schrodinger-scenario 
mately related to the reciprocity relations of the preceding Section, as discussed in [30| and 
states that 

The potential function U{x) belongs to £2 and thus admit of a series expansion in 
X, x^, x^,etc. 29|. The Ak{x) themselves belong to £2 as well and can be series-expanded 
in similar fashion. This enables us to base our future considerations on the assumption that 
the a priori knowledge refers to moments of the independent variable, i.e., 

(A) = {x') , (14) 

and that one possesses information on M moment-mean values {x^). Our "information" 
potential U then reads 

U{x) = -W>^k^'- (15) 
^ k 

and Eq. (|T3|) allows one to immediately obtain 



(^) = 8 E^A, {A,)- {A,=x'), (16) 

and thus, via (fT6l) and the above mentioned relation I = — 4:{dxx), a useful, virial-related 
expression for Fisher's information measure can be arrived at. 



/ = - E^^fc(^')' (17) 
k=l ^ 

which is an exphcit function of the M physical parameters (x^) and their respective Lagrange 
multiphers A^. Eq. (fT7|) encodes the information provided by the Virial theorem. Thus, we 
have two different ways of expressing /, namely, ([8]) and f|T7|) . Interesting things happen 
if we put them together. Since is given by (|TTi) as [dl /d{x^)], inserting the reciprocity 
relations ( fTTl) into ( |T7l) we are led to 



Eq. (fT8|) constitutes an important result, since we have now at our disposal a differential 
FIM-equation. Dealing with it should allow us to find / in terms of the {x^) without passing 
first through a Schrddinger equation first, a commendable achievement. This is a linear par- 
tial differential equation that an extremal / must necessarily comply with. This constitutes 
one of the main present results. It is not clear (yet) whether from such an J— form we can 
extract an amplitude ip satisfying a Schrodinger equation. Our I could however be related 
to an approximate solution to Schrodinger's equation. For convenience we now recast our 
key relations using dimensionless magnitudes 

where [/] and [{x^)] denote the dimension of I and {x'^), respectively. Thus, the differential 
equation that governs the FIM-behavior, i.e., (1T8|) . can be translated into 

^ k BT 

^=-T.7^(^k) TT7T7T , I = I{m,---, (A-m)), (20) 

fc=l ^ o{^k) 

which is a first order linear nonhomogeneous equation with M independent variables. All 
first order, linear partial differential equations (PDEs) possess a solution that depends on an 
arbitrary function, called the general solution of the PDE. In many phy sical situations this 
solution if less important than other solutions called complete ones |32|-|34| . Such complete 
solutions are particular PDE solutions containing as many arbitrary constants as intervening 
independent variables. Let us look now for a special complete solution of our PDE (!20ll . 
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whose usefulness will be illustrated via two physical examples below (the treatment of general 
solutions is postponed to Section VI). We first set 



M M 

1= ^ Xfc = exp[^((A'fc)) 

fc=l k=l 



and substituting (ETil into (ITHi) leads to 



k=l 



The above relation entails 



(21) 



(22) 



9iW) = -Y InlWI+Cfc 



k {Xk) 

where Ck is an integration constant. Finally, substituting ( 123|) into (12T]) we arrive at 



(23) 



M 



^=E C'fc exp(^- ^ InKA-fc)!) , 



C.= e"^- > 



which can be recast as 



M 



X((A'i),...,(A'm))= E C'fc l(A'fc) 



-2/fe 



/c=l 



or, in function of the original input-quantities f fTOjl 



M 



/((xi),...,(x^))= 5: (a:'=) 



-2/k 



k=l 



an intriguing result. We enumerate below the main properties of this minimal I. 



(24) 



(25) 



(26) 



• FIM-domain 

Obviously, it is 



Dom[I] = {{{x'),...,{x''))/{x') e ^o} 



FIM-monotonicity 

Differentiating ( l26l) we obtain 
dl 



Ik 



k {x'^) 



-2/k 



(27) 



Therefore, if {x^) > , / is a monotonically decreasing function in the (x'^) -direction. 
Also, for (x'^) > , from the reciprocity relations f|TT]) we have, 

Xk = ~ ^Ck{x')-^'+'^/' < 0. (28) 

FIM-convexity 

This is a necessary property, since the entropy is concave. By differentiation of the 
expression (1271) one obtains 

^ , 4 -2{l + k)/k 



d{x'')d{x^) \ 2/ F 
from which we can assert that the Fisher measure is a convex function. It is then 
guaranteed that the inverse of dkdjd exists. 



IV. THE REFERENCE QUANTITIES Ck 

FIM is an estimation measure known to obey the Cramer Rao-bound ([2]) j4|. The best 
estimator exhibits a CR relation as close to unity as possible. Thus, the reference quantities 
Ck should be chosen in a manner that respects this condition. Here we are interested in 
simple situations that illustrate the concomitant procedure. More involved situations will 
be treated elsewhere. 

Since the reference quantities Ck contain important information concerning the reference sys- 
tem with respect the which prior conditions are experimentally determined, it is convenient 
to start by choosing an appropriate reference one. 



Minimum of the information potential 

We consider it reasonable to incorporate at the outset, within the /—form, information 
concerning the minimum of the information potential U{x). Assume that this information- 
potential 

U{x) = - -Y.>^kx\ 

^ x=l 

achieves its absolute minimum at the "critical point" a; = ^, 

u\0 = , U^in = f/(0- (30) 
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Effecting the FIM-translational transform u = x — ^ leads us to 

/ = 

with (see the Appendix) 



k=l ^ 



k=i ^ 



(31) 



A! 



{uy = {{x - 0') 



(32) 



where U^''\^) is the fc*^ derivative of U(x) evaluated at a; = ^ and ( )' indicates that the 
pertinent mean value (x— moment) is evaluated for translation-transformed eigenfunctions. 
The corresponding FIM-explicit functional expression is built up with the A^— non-vanishing 
momenta (A^ < M) {{u^)' ^ 0) and is given by 



N 



k=2 



1= T.Ck {u'y 



-2/k 



N 



k=2 



E Ck {{x - 0") 



-2/k 



(33) 



where we kept in mind that A* = —8U'{^) = 0. A glance at the above expression suggests 
that we re-arrange things in the fashion 



/ = C, {{x - 0') 
Taking now into account that 



N 



E Ck {{x - 0') 

k=3 



-2/k 



(34) 



{x-O = 

{{x - 0') = {x') - 2ax) + e 



we get 



a 



N 



{{x - 0') = {x') - {xf = a' 



(35) 



+ Y.Ck {{x-o") 



from which we obtain 



fc=3 



N 



-2/k 



(36) 



Ia'= C2 + a^E Ck {{x-O") 

k=3 



-2/k 



> 1 



(37) 



Therefore, / preserves the well-known Cramer- Rao /—bound js^ / o"^ > 1. The above seems 
to indicate that if no moment > 3 is a priori known, the lower bound can be reached for 
C2 = I. For A; > 3 additional considerations apply that will be discussed elsewhere. 
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V. TWO PHYSICAL EXAMPLES 

So as to illustrate the above considerations we are going to consider two simple and instruc- 
tive examples. We take the mass m = 1 and h = 1. 

Harmonic oscillator (HO) 

The prior information is given by 

{x^} = ^ , M = 1 , k = 2. (38) 

The minimum of the potential function obtains at the origin ^ = 0, 

The pertinent FIM can be obtained using f l33p with u = x — ^ = x, 

I = Ii{x')) = C2 {xy , 

and, the CR bound is saturated when C2 = 1, 

/ = C2 = 1 =^ I = {xY' . (39) 

The corresponding Lagrange multiplier can be obtained by recourse to the reciprocity rela- 
tions (fni) and (1391), 

a|) — <^^>" ■ («' 

The prior-knowledge (1551) is encoded into the FIM (15^ . and the Lagrange multiplier A2 (HO]) . 

/ = -1 = 2u;; A2 = - {x^) = - Aoo^ . (41) 

and the a— value can be obtained from ( |TOl) . 

a = I - X2 (x^) = 4 (42) 

as we expect. 
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Harmonic oscillator in a uniform field 



We consider a charged unit-mass particle moving in the HO potentiaL The electrical charge 
is q and there is a uniform electric field e, in the a;— direction. Our prior knowledge is given 



by (29| 



We look first for the ^-point at which U{x) is minimal. 

U{x) = — (Ai X + X2 X 



2 



f/'(0 = -^(Ai + 2A20 = ^ = -Yk ■ ^^^^ 

The translational transform u = x — ^ implies that 

w = {x-o = {^)-^, i^y = {{x - 0') = ix') - 2^{x) + e . (45) 

The translation transformed FIM is now given by 

I = C2 {u')'-\ (46) 
and, the CR bound is saturated when C2 = 1, 

/ {uy = C2 = 1 =^ I = {u^-' . (47) 
The reciprocity relations lead us to 

o{x) d{u^)' d{x) 

^ ^ Wyy 'di^ = - ) • (49) 

From the prior knowledge ( 143|1 and using (HSll we have 

(x) = e = ^ , (50) 

LiJ 
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then, inserting ([50]) and (jH]) into dH]) - we get 



(52) 



X2 = -{U')'-' = 



2 11 (2.)^ 
- (2 c.)^ 



8 q e 



(53) 
(54) 



The corresponding translational transform a— value can be obtained substituting ( l52i) -(l& 
into (Uni), 



a=I - X,{x) - X2{x') = 2u-SqM + Au^ ( ^ + C^y"] = 4u; , 



(55) 



and the corresponding a— value is given by (see Appendix), 



a = a + 8 f/(0 = 4w - 4 



(56) 



as we expect. 



VI. GENERAL SOLUTION OF THE DIFFERENTIAL FIM-EQUATION 



We discuss here this issue for the sake of completeness. Our F/M— equation is a first 
order linear nonhomogeneous differential equation. We are following 32|-|34| in looking for 
the general solution. For a first-order PDE, the method of characteristics allows one to 
encounter useful curves (called characteristic curves or just characteristics) along which the 
PDE becomes an ordinary differential equation (ODE). Once the ODE is found, it can be 
solved along the characteristic curves and transformed into a solution for the original PDE. 
The characteristic system of Eq. (l20l) is 



leads (for {X^) ^ 0) to 



dl 
T 



(57) 



d{Xk) 



d{X,) 



{k/2){Xk) {1/2){X,) 



- \n\{Xk)\+Ck 



2 In 1(^-1)1+01 



In 



2/k 



,2/k 



In 
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dX d{Xi) 



X (l/2)(A'i) 



hk_,^e^--^^ = \{X,)Wx^)r^'^ (58) 

ln|X| +cx = -2 ln|(;fi)| + Ci 

In [e'^ |X|] = In [e"^ \{Xi)\-^' 

e"^ \X\ = e"^ \{X^)r^ 

i 

bM = e''-'^ = \{X,)f \X\ . (59) 



We have now constructed an integral basis for the characteristic system (j5j 

h = u,{{X,), {Xm),X) , . .. , bM = um{{Xi), {Xm),X) , (60) 
and the general solution of equation (1201) defined as 

<I>(mi,M2, . . . ,um) = 0, (61) 

is given by 

<i.(i(A'i)n(A'2)r\---,i(A'i)n(A',)r^/'^ ,...,\{x,)f\{XM)\-'/'' \x\) = o, 



where $ is an arbitrary function of the M variables. Solving this equation for X yields a 
solution of the explicit form 

where is an arbitrary function of (M — 1) variables. 

Cauchy problem and the existence and uniqueness of the solution to our PDE 

One of the fundamental aspects so as to have a useful PDE for modeling physical systems 
revolves around the existence and uniqueness of the solutions to the Cauchy problem. Here 
we show that such requirements are satisfied by our pertinent solutions. We start by casting 
(l20l) in the normal form 



^ f({Xi),---AXm),X,J^,...,-^] (63) 
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where 



Fx = F[{X^),■■■,{XM),Z, 



dX 



dX 



d{X2y 'd{XM) 



^ k dX 



(64) 



and we see that F is a real function of class in a neighborhood of 



dX 



d{Xu) 



4-1, A; = 2,---,M 



Then, if ■?/'((A:'2), {Xm)) is also a function of class such that 

dil) 



^(^l,--,^Af-l) = c 



d{Xk 



4-1, A; = 2,...,M. 



(65) 



exists a solution X of (!63|) in a neighborhood of {Xi) = a and {Xk) = ^k-i that satisfies 



X(a, (A'2),---,(A'm))=^((A'2),---,(A'a/)) 



(66) 



and is of class C^. 



Regarding Cauchy-uniqueness, it is known that if F satisfies the Lipschitz condition [35|, 



M 



k=2 



dx' dx 



+ K2 



d{Xk) d{Xk) 
F^, =f({X,),---,{Xm),X 



X -X 

, dX' 



Ki, K2 = const. (67) 
dX' \ 



d{X,)' 'd{XM)J 

then, the solution of the initial value problem for f l63|) is unique. Note that in our case the 
above condition is verified always since the Legendre structure the theory guarantee that 



dX dl 

oc 



d{Xi) d{x) 



Ai < 00. 



(68) 



VII. CONCLUSIONS 

It was known for some time js] that, for one-dimensional scenarios, a minimal Fisher's infor- 
mation measure I is associated to a Schrodinger's equation in which the role of the potential 
is played by a weighted sum of the a priori known expectation values of M quantities (x*^). 
The weights are the corresponding Lagrange multipliers A^. Legendre reciprocity relations 
linking J, the A^ and the (x*^) apply. 
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We have here demonstrated that insertion of virial theorem-tenets into this Legendre struc- 
ture leads to a differential equation for /. The equation is analytically solvable and its 
solution provides us with explicit new expressions for I in terms of the input-information 
contained in the M expectation values {x^). In other words, we can directly codify the in- 
formation provided by such set of expectation values in an J— form without previous appeal 
to a Schrodinger equation. Additionally, this (partial) differential equation (PDE) may be 
viewed as a new constraint that the solutions ipn of Schrodinger's equation (SE) must satisfy. 
In this regard, our PDE may lead to useful criteria for checking approximate treatments of 
SE. If (f)n is an approximate SE-solution, the fact that the functional /[0„] verifies the PDE 
would constitute an indicator of the quality of the approximate solution 0„. 

An application to simple examples has illustrated these considerations. Of course, as is the 
case in the MaxEnt environment, the usefulness of fl26l) depends on how adequate is our 
input information for describing the situation at hand. Maximal entropy or minimum FIM 
are just the best ways to exploit that knowledge. 



Appendix: Translation transform of FIM 

The potential function 



1 M 



U{x) 



can be Taylor-expanded about x = ^ 



M 



u{x) = y: 



k ! 



{x-0' 



k=0 



The translational transform u = x — ^ leads to 



M 




u{u)= u{u + o= E 



(69) 



fc=0 



which can be recast as 



Uiu) = - - E Ku' , 



(70) 
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with 

= - 8^^ = -|^f: j(j-l)(j-2)---(j-A: + l)A,e-^ (71) 

■ i=i 

The FIM-translational transform u = x — C, is obtained from ([HD in the fashion (( )' indicates 
that the moment is calculated for translation-transform eigenfunctions) 



/ . _ 4 / = - = - 4 ) , (72) 



where %() = ip{u) is the translation transform of il){x). Now, using the translation transform 
of ([6]) one easily finds 

M 



I ^nia + Y. \l u''\ ipn du , (73) 

V fc=0 / 



and one realizes that 



M M 

I = a + Y: K{u'y = « + E Kiu'Y , (74) 

k=0 k=l 

where 

a = a + Ao = a-8f/(0. (75) 
Also, the virial theorem ( lT3l) leads to 



The translation-transformed moments (m'^)' are related to the original moments as 



{uy = j ^\u) du = j + i)du = j{x- 0'' ^'(a;) dx = {{x - O'") 

By recourse to the Newton-binomial we write 

J{x- O'^ ^'(x) dx = J2 (-1)' ( ^ j / ^'(^) dx, (77) 
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and then we finally have 



3 



(78) 
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